A. Thiaville\ Y. Nakatani^'^, J. Miltat^ and Y Suzuki^ 

^ CNRS-Universite Paris-sud, 
Laboratoire de physique des solides, 
Bat. 510, 91405 Orsay cedex, France, 
^ University of Electro- communications, 
Chofu, 182-8585 Tokyo, Japan, 
^ Osaka University, Dep. of Materials Engineering, 
Toyonaka, 560-8531 Osaka, .Japan 

(Dated: January 21, 2005) 

In order to explain recent experiments reporting a motion ol magnetic domain walls (DW) in 
nanowires carrying a current, we propose a modification of the spin transfer torque term in the 
Landau-Lifchitz-Gilbert equation. We show that it explains, with reasonable parameters, the mea- 
sured DW velocities as well as the variation of DW propagation field under current. We also 
introduce coercivity by considering rough wires. This leads to a finite DW propagation field and fi- 
nite threshold current for DW propagation, hence we conclude that threshold currents are extrinsic. 
Some possible models that support this new term are discussed. 

PACS numbers: 72.25.Ba, 75.70.Kw, 72.25.Pn 



Recent research on magnetic nanostructures has shown 
that effects caused by an electric current flowing across 
a nanostructure may dominate over the effects due to 
the field generated by the same current P, Ql • Most of 
the work up to now, experimental and theoretical, has 
been devoted to the 3-layer geometry (2 magnetic lay- 
ers separated by a normal metal spacer, with lateral di- 
mensions well below the micrometer so as to get single 
domain behaviour). Under current, generation of spin 
waves layer switching Q and precession of magne- 
tization y have been observed. All these results could 
be qualitatively explained by the spin transfer model I'l . 
On the other hand, the situation with an infinite number 
of layers, namely a magnetic nanowire containing a mag- 
netic domain wall (DW), has just started to be studied. 
Here, under the sole action of a current, the DW may 
be moved alongthe wire, as confirmed by several exper- 
iments 0, 0, US) IM Q • The situation is however more 
complex than with 3 layers, as the evolution of the cur- 
rent spin polarization across the DW has to be described, 
so that a deep connection with the problem of DW mag- 
netoresistance exists. Two limits have been identified by 
the theories, namely the thin and thick DW cases. The 
length to which the DW width has to be compared is, 
depending on the model, the s pin diffusion length 
the Larmor precession length |l2| . or the Fermi wave- 
length |0| . These lengths are below or of the order of a 
nanometer in usual ferromagnetic metals. In the exper- 
iments cited above (nanowires of typical width 100 nm 
and thickness 10 nm), the DW width is of the order of 
100 nm, much above these lengths. Thus, the approxima- 
tions of current polarization adiabaticity and full transfer 
of angular momentum to the local magnetization apply, 
and the Landau-Lifchitz-Gilbert (LLG) equation of mag- 
netization dynamics becomes 

TO = 7o-ff X TO -f- am xm— (m-VIto. (1) 



The symbols are: to unit vector along the local mag- 
netization, 70 gyromagnetic constant, H micromagnetic 
effective field, Ms spontaneous magnetization and a 
Gilbert damping constant. The velocity u is a vector 
directed along the direction of electron motion, with an 
amplitude 

u = JPgfiB/ {2eMs) , (2) 

where J is the current density and P its polarization 
rate {u is positive for P > i.e. for carriers polar- 
ized along the majority spin direction). For permalloy 
the factor 5^B/(2eA/s) amounts to 7 x 10"" m^/C. The 
local form (HI, already introduced by several authors 
19, ,20, 21, 22], can be seen as the continuous limit of 
the Slonczewski spin transfer term between infinitely thin 
successive cross-sections of the nanowire. 

The consequences of ^ on the dynamics of a trans- 
verse wall (TW - Fig. ^) have been well studied quali- 
tatively [13, Il4|. Using the integrated equations for the 
DW dvnamics |l5l[l6']. it was found that 2 regimes exist. 
At low current, the spin transfer torque is balanced by 
an internal restoring torque. The DW magnetization tilts 
out of the easy plane (as first explained by Berger ^1^), 
but the DW does not move steadily under current. Above 
a threshold, the internal torque is not sufficient and DW 
motion occurs, together with a continuous precession of 
the DW magnetization (in Fig. ^ a d ^ a d..., the 
bar meaning tt rotation around the x axis). The thresh- 
old is linked to the steepness of the potential well that 
defines the equilibrium orientation of the TW magneti- 
zation (mainly a magnetostatic effect, see Fig. 

These qualitative findings were fully confirmed by mi- 
cromagnetic computations [l^ uncovering, however, a 
big discrepancy with experimental results. The intrin- 
sic current thresholds computed for steady propagation 
of a TW in a perfect wire were found to be more than 
10 times larger than experimental values. Thermal flue- 
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FIG. 1; Schematic domain wall structures in a soft nanowire. 
They are governed by magnetostatics and depend on the wire 
cross-section: (a) transverse wall, (b) perpendicular trans- 
verse wall, (c) vortex wall. As the fabricated wires are very 
wide compared to the exchange length (5 nm for permalloy), 
the energy maximum situation (b) decays to a wall with an 
antivortex (d) |24|. In cases (a, b) the DW magnetization 
angle <f> with respect to the y axis is and 90°, respectively. 
The stable structures are (a) at smaller and (c) at larger cross- 
sections |2a|. 



tuations cannot change the situation because the energy 
barrier against DW magnetization rotation proves very 
large (for a 120 x 5 nm^ permalloy wire, the energy differ- 
ence between states (d) and (a) in Fig.n]is 230 kBTamb)- 
Therefore, we need to introduce corrections to perfect 
adiabaticity and pure local spin transfer in the micro- 
magnetic models. In this letter, we show how the inclu- 
sion of a second term describing the current torque, as 
introduced previously in the 3 layers geometry 2311, can 
resolve this discrepancy. Our approach is twofold: first 
we study how micromagnetic results are modified when 
this term is phenomenologically included. In a second 
step we examine several models of the current effect that 
may justify the existence of such a term, and try to esti- 
mate its magnitude. 

The modified LLG equation, again obtained in the con- 
tinuous limit, now reads 



m 



(3) 

As discussed later, the non-dimensional parameter /3 is 
much smaller than unity, hence comparable to a. The 
solved form of ^ shows that the new term modifies the 
initial spin transfer torque by a second order quantity, 
and directly competes with the damping term associated 
to the first term. This is very appealing, as one sees 
directly from Q| that for zero damping the solution under 
current is just the solution at zero current, but moving 
at velocity u. 

In order to test the influence of this new term, micro- 
magnetic numerical simulations were performed. Simi- 
larly to previous work a moving calculation region 
centered on the DW represented an infinite wire. Wire 
width and thickness corresponded to experiments 0, ll| , 
and material parameters were those typical of permalloy 
(Ms = 8 X 10^ A/m, exchange constant A = 10"^^ J/m 
and no anisotropy) with, in most cases, damping fixed to 
a = 0.02 The mesh size was 4x4x5 nm'^. It was 
checked first that the Oersted field generated by the cur- 
rent had a negligible effect for current values correspond- 



ing to experiments. Indeed, because of the small thick- 
ness of the sample, this field is essentially perpendicular 
and localized at sample edges. Its maximum value for an 
extreme case u = 100 m/s, P = 0.4 in a 120 x 5 nm^ wire 
is ^jLqHz = 0.017 T, a value well below the perpendicular 
demagnetizing field (1 T), that results in a very small 
out of plane magnetization rotation. Consequently, all 
results shown below were computed without the Oersted 
field. Transverse, but also vortex walls, were considered. 

Figure shows the wall velocity v under zero field as 
a function of the velocity u (proportional to the current) 
with /? as a parameter. The sample is 120 nm wide and 
5 nm thick (2|, with a DW of the transverse type (Fig^). 
The /3 = curve displays an absence of DW motion for 
u < Uc — 600 m/s, and a rapid increase of v towards u 
above the threshold, as inferred previously [l3. l2^. 
Curves with f3 > differ markedly from this behaviour, 
as DW motion is obtained at any finite u (keep in mind 
the wire is perfect). Velocity increases linearly with u and 
f3, up to a breakpoint where it decreases (or increases 
if /3 < a - not shown). Results obtained for other a 
values indicate that all of them can be described as v = 
(3u/a below the breakpoint. Above breakpoint, the DW 
structure examination reveals the periodic injection of 
antivortices at wire edges, which cross the wire width 
and are expelled, similarly to the field driven case or 
to the current driven case with /? = above their 
respective thresholds. 

The same behavior is obtained at low u for a VW struc- 
ture in a larger wire (width 240 nm, thickness 10 nm). 
This surprising result (for field-driven propagation, wall 
mobility is much lower for a VW j2Sj) can be explained 
by the global arguments of Thiele [2y| about steady-state 
motions. Indeed, Eq.|21can be transformed to a general- 
ized Thiele relation 



F + Gx{v-u)+ D {av- l3u)^Q, 



(4) 



with F the static force, G the gyrovector (along z here) 
and D the dissipation dyadic |2^. The solution is then 
V = (3u/ a as long as the restoring force keeping the vortex 
in the wire can balance the gyrotropic term. Otherwise 
the VW transforms into a TW by lateral expulsion of the 
vortex. 

Experiments on single-layer wires of permalloy 0, 0, 
llo| show however that a threshold current exists for DW 
propagation, similarly to field-driven propagation (for ex- 
ample, a propagation field Hp ~ 20 Oe was measured in 
Ref. 13). This is not reproduced by our results for per- 
fect wires. A roughness of the wire edges was therefore 
introduced, characterized by a mean grain size D (see 
Ref. [2^ for details). Fig|5|3 shows how the results of 
FigEt are modified for D — IQ nm, a value which re- 
sults in an extrapolated propagation field H-p ^ 10 Oe, 
with velocity vs. field curves given in Ref. [ig. At low 
(3 (< 0.02) the DW is blocked for u < Uc, showing that 
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FIG. 2: Steady velocity computed for a transverse domain 
wall by micromagnetics in a 120 x 5 nm^ wire as a function of 
the velocity u representing the spin polarized current density 
@, with the relative weight (/9) of the exchange field term as 
a parameter. Open symbols denote vortices nucleation. The 
shaded area indicates the available experimental range for u. 
(a) perfect wire and (b) wire with rough edges (mean grain 
size D — 10 nm). Dashed lines display a fitted linear relation 
with a 25 m/s offset. 



roughness-induced DW pinning dominates. Above Uc, 
the DW moves via nucleation of antivortices and vor- 
tices, again similarly to the perfect case with (3 — 0. If 
(3 is larger (0.02 — 0.1), a window of DW propagation 
without vortex nucleation exists (filled symbols) , includ- 
ing values of u well below Uc- Larger /? values result in 
DW motion at lower u. However, with our coercivity 
model and the absence of thermal fluctuations, continu- 
ous DW propagation at "low" velocities {v < 300 m/s) 
is not possible as the DW eventually stops at positions 
with a larger wire width fluctuation. This may be largely 
due to our coercivity model, as in the same conditions a 
VW moves already at w « u = 100 m/s, and for regu- 
lar notches TW propagation at w < 10 m/s is obtained. 
By comparing with Fig. it is clear that the rough 
wire data cannot be fitted by merely introducing an off- 
set Up into the results for a perfect wire, as done usually 
in the field-driven case. On the other hand, simple pro- 
portionality of w to M gives a very good fit. However, it 
would be as good if a small offset existed, as shown by 
the dashed lines in Fig. ^jp that assume up = 25 m/s . 



Thus, this model of coercivity does create a propagation 
field Hp and is consistent with a threshold current up for 
DW propagation, even if the data do not allow to deter- 
mine precisely its value. An important additional result 
of the simulations is that, for a rough wire, the velocity 
is reduced when compared to the perfect case relation 

V = ((3/a)u. The reduction factor is not a constant, as 
for example at Z? = 10 nm, a = 0.02 and f3 — 0.04, the 
factors for a TW are 0.6, 0.7 and 0.4 for 120x5, 240x5 
and 240x10 nm^ strips, respectively, whereas they are in 
the range 0.5 — 0.6 for a VW. 

Experimental propagation threshold currents up 0, H, 
{u\ for permalloy strips amount to P x 50 m/s (TW) 
and P X 80 m/s (VW), for cross-sections 120x5, 240x10 
and 200x(5-35) nm^, the current polarization P being un- 
known. These values are compatible with our results at 
_D = 10 nm, which was chosen in order to reproduce typi- 
cal propagation fields. Although the value of (3 cannot be 
determined at this stage, we infer from our results that 
it should be larger than 0.01. In addition, the measure- 
ment of the propagation field under current showed 
a linear variation of Hp vs. u, with AHp = 3 Oe for 
Au — P X 100 m/s. As for a perfect wire [l^ one has 

V = fiH + (3u/a, we evaluate AHp/ Au = (3/{aiJ,). The 
mobility calculated by micromagnetics for this nanowire 
being /i = 26.5 m s~^Oe~^ H^, we estimate, for a = 0.02 
and P = 0.4, (3 = 0.04. A micromagnetic calculation for 
this situation, on the same 120x5 nm^ wire with a — 0.02, 
/3 = 0.04 and u = ±50 m/s gives AHp 7.8 Oe which 
fully agrees with the above (0.4 x 7.8 = 3.1). Another ex- 
perimental data is the DW propagation velocity v (only 
one measurement, Ref. 0). It was seen to increase from 
2 to 6 m/s as the current velocity was increased from 
P X 80 to P X 90 m/s. The average slope of these data 
is consistent with, in a perfect wire, /3 = 0.03, a = 0.02, 
and P — 0.4. Note that (i) the (3 value should be larger 
because of the velocity reduction in a rough wire, and 
(ii) the determination of /3 is subject to uncertainty, as 
the coercivity source we considered may not be the only 
meaningful one. To summarize, we see that the exper- 
imental results can be reproduced by introducing the /3 
term in JSJ, with a value of about 0.04 for permalloy. 

We describe now how the numerical micromagnetic re- 
sults shown above can be qualitatively understood in the 
TW case, in the framework of the extended Slonczewski 
equations describing DW motion. These are constructed 
by assuming a given DW profile with some parameters 
(here the wall position q and magnetization angle - see 
Fig^. In the case of pure spin transfer the LLG equation 
could be derived from a least acti on p rinciple, allowing 
an easy derivation of the equations [13 . This is no longer 
true with the additional term and the original procedure 
must be applied 0, 0| . For a one-dimensional TW pro- 
file, characterized by a DW width parameter A, one then 
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FIG. 3: Domain wall propagation velocity under current as 
obtained from the one-dimensional model, assuming a perfect 
wire and a constant wall width parameter A, for a damping 
constant a = 0.02. The curves' shape is extremely similar to 
the numerical results (Fig. I^Ji)- Velocities are normalized to 
the Walker threshold ww = 7oAii"K/2. 



the latter, dominant for thin walls, enters the integrated 
DW dynamics equations © similarly to our second term. 
Thus, for a DW of finite thickness the second term should 
remain present in some proportion. As the characteris- 
tic length is here the Fermi wavelength |l3l |. we expect 
/3 <C 1. Note however the lack of an expression in the in- 
termediate DW thickness regime that would predict the 
value of /?. 

Zhang et al. calculated spin accumulation for a current 
flowing across a multilayer, and the resultant torque on 
the magnetization, with the inclusion of the s-d exchange 
interaction [ill Eof . Extending their model j^J to the 
case of a linearly varying magnetization gives both terms 
of with a slightly modified equivalent velocity u and 
a P factor given by 

(3 = {Xj/Xsff = h/ (Jnj) , (7) 



gets 



7oi?a+/3u/A ; q/A-aq 



with J the s-d exchange interaction energy, Tgf the spin- 
7oi?K sin0cos0-|-u/zfl,ip time, Aj and As/ the associated diffusion lengths. 



(5) 

with A = A(0) the wall width parameter mini- 
mizing the DW energy at wall magnetization angle 
(j). In the simplest ID model one has A((/)) = 

1 /2 

(A/ (Kq + Ksin^ (t))) , with Kq the effective longi- 
tudinal anistropy and K the transverse one (i?K = 
2K/ (fioMs)). One sees from (0) that the additional 
term enters the equations similarly to an applied field 
Hcquiv — ^P/iloA), explaining DW motion at any u in 
a perfect wire. At zero applied field one gets 

• a \ P — a u JqHk 



1 



■ sin 2(j) 



q = (3u/a—A(f)/c 



a A 2 

For |w| < (7oi?K/2) A(7r/4)a/|/3 - a| |2g a stationary 
regime exists with q = j3u/ a. The DW velocity is remark- 
ably independent on the DW width, as already shown 
with Q). Above the threshold oscillates and so does 
the velocity. The average wall velocity, plotted in Figl^l 
is very similar in behaviour to numerical results (Fig|2^) . 
The non-trivial effect of roughness is not included in this 
ID model, but if one takes it into account by an increased 
dissipation (larger a) a velocity reduction in rough wires 
is natural. 

We finally turn to discussing possible origins of EqOl 
One may first remark that the new term exhausts the 
possibilities of introducing linearly the magnetization 
gradient along the electric current direction into the LLG 
equation. 

The approach of Tatara and Kohno put forward 
two effects of an electric current across a DW, called spin 
transfer and momentum transfer. The former, dominant 
for thick DWs, has the same form as our first term while 



Note that with Aj = 1 nm and Xgf — 5 nm for permal- 
loy, one gets = 0.04. As for pure metals As/ can be 
much larger [2(j| , one expects that (3 may become smaller 
and, hence, that DWs may prove more difficult to move. 
This should be investigated by experiments. One should 
also be aware that the time Ts/ in ((TJ may not be ex- 
actly the spin-flip time. This time appears in fact as the 
transverse relaxation time in the evolution equation of 
the spin accumulation. If part of this relaxation occurs 
through interaction with the localized moments, for ex- 
ample by incoherent precession around the local magne- 
tization, then the moment is transferred to the localized 
spin system and does not contribute to /?. 

In conclusion, we have shown that the inclusion of 
a second term describing the spin transfer torque, akin 
to the exchange field term considered previously, might 
solve the puzzle of the quantitative explanation by micro- 
magnetics of the current induced DW motion 31]. This 
view is supported by a direct comparison of experimen- 
tal results (DW velocity, reduction of DW propagation 
field under current) to micromagnetic simulations of re- 
alistic nanowires (excepting pinning phenomena). The 
model predicts that DWs should move under smaller cur- 
rents for more perfect wires (having a lower propagation 
field). As our approach was phenomenological, a basic 
understanding of the form of the spin-transfer torque in 
continuous structures as wide as a domain wall is still 
very much needed. The magnitude of the new term, and 
its material dependence, should allow to discriminate be- 
tween models. 

A.T. acknowledges fruitful discussions with F. Piechon 
and N. Vernier. 
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